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Abstract. I show that the spectrum and morphology of a recent Fermi-LAT observation of 
the Galaxy center are consistent with a millisecond pulsar population in the nuclear Central 
stellar cluster of the Milky Way. The Galaxy Center gamma-ray spectrum is consistent with 
the spectrum of four of eight globular clusters that have been detected in the gamma-ray. A 
dark matter annihilation interpretation cannot be ruled out, though no unique features exist 
that would require this conclusion. 
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1 Introduction 

The ability of gamma-ray observatories to provide a window on dark matter annihilation 
signals has been known for some time, (e.g. [1]). The recent manuscript by Hooper & Goode- 
nough [2] performs a detailed analysis to extract the signal in the gamma-ray of the dynamical 
center of our Milky Way Galaxy as observed by the Large Area Telescope (LAT) aboard the 
Fermi Gamma-Ray Space Telescope. That work interpreted the morphology and spectral 
feature of the signal within the inner 1.25° (175 parsec radius) from the Galactic dynamic 
center to be inconsistent with any known astrophysical sources, and further claimed for it to 
be accounted for by the presence of annihilating light (~8 GeV) dark matter. In this short 
note, I show that the signal as seen by Hooper & Goodenough (hereafter HG) is from the 
Galactic Central stellar cluster and is consistent with prior observations of gamma-ray emis- 
sion from massive stellar globular clusters with a population of millisecond pulsars (MSP) 
[3]. 
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Figure 1. Shown is the spectrum E 2 dN/dE of the HG source (thick black points with errors), and 
its best fit power law with exponential cutoff curve (blue). The HG spectrum power-law index is 
consistent with the globular cluster gamma ray emission of Omega Cen (green), NGC 6388 (cyan) 
and M 28 (magenta) [3]. The spectrum of the Geminga pulsar (light grey) [4] has a nearly identical 
peak energy and exponential cutoff and a softer power-law index than the HG source. The spectra 
of Omega Cen and M 28 are shifted by +0.14 dex and +0.34 dex in energy, respectively, to illustrate 
the consistency in power-law index. NGC 6388 is not shifted in energy. The HG source is consistent, 
within errors, with the shape of the intrinsic spectrum of all the plotted globular cluster sources. All 
sources are normalized to the HG peak flux. 

2 Spectrum 

In general, pulsar gamma-ray spectra have a flux peak at ~ 2 — 3 GeV and a spectrum 
dN/dE oc E~ T exp(E/E c ) with a power law with index V in the range 0.5-2.5 and exponential 
cutoff energy E c at 1-5 GeV, a range consistent with the HG source [5]. Note that HG compare 
to the average millisecond pulsar spectrum and find it is not consistent, and then reach the 
conclusion that no astrophysical source can produce the HG spectrum. However, this is not 
appropriate, but rather consistency with any set of pulsars or stellar cluster population of 
pulsars must be ruled out before claiming inconsistency with astrophysical sources. 

The spectrum of the HG source can be well fit by a pulsar-like power law plus exponential 
cutoff spectrum. The spectrum below 7 GeV is well fit by dN/dE oc E~ r enp(E/E c ), with 
T = 0.29 ±0.12 GeV and E c = 1.34 ±0.12 GeV and a reduced r7 D OF = 0.5. The spectrum 
and best fit curve is shown in Fig. (1). 

Several pulsars in the First Fermi-LAT Catalog of Gamma-ray Pulsars [5], including 
J1958±2846, J2032±4127 and J2043±2740, have a power-law index and exponential cutoff 
consistent with the HG source. The Geminga pulsar has a detailed published gamma-ray 
spectrum which has a similar spectral peak, 2 — 3 GeV and cutoff as the HG source [4] . The 
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Geminga spectrum is shown for comparison in Fig. (1). However, Geminga has a shallower 
power-law spectral index T = 1.08 ± 0.02. 

However, it is not expected that the extreme density of the Galactic Central stellar 
cluster harbors only a single pulsar. The Central stellar cluster (nuclear bulge) of the Milky 
Way is coincident with the observational region of HG central spectral signal, and extends 
across the full ~ 1° from the dynamical center (or ~300 pc across) of that spectral analysis 
region. The spectrum of the HG Central stellar cluster region (r = 0.29 ± 0.12, E c = 
1.34 ± 0.12) is consistent with four of the eight globular clusters that have been observed by 
Fermi-LAT and seen to have gamma-ray emission [3]: 
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Three of these spectra are shown in Fig. (1). Note that a direct quantitative comparison of 
the fit to the HG source spectrum presented above with the globular cluster fits, Eqs. (2.1), 
cannot be rigorously performed since the sets of fits use different methods. The fit performed 
in this analysis is to the points and errors presented in HG, while the fits in Eqs. (2.1) from 
Ref. [3] are performed with the Fermi-LAT analysis tool gtlike, which bins the spectrum 
differently and more accurately incorporates the photon statistics. The fit between HG and 
the globular clusters are obviously qualitatively consistent (Fig. 1), but a detailed statistical 
comparison can only be established with gtlike analyses of all regions. Most importantly, 
over-subtraction or under-subtraction of the diffuse component of an observation can lead 
to a systematic uncertainty that must also be included in a full error analysis of the Galaxy 
Central stellar cluster spectrum [6]. This was not performed in the HG analysis. Over- 
subtraction of the diffuse component may lead to a high estimate of the intrinsic spectral 
index T. 



3 Morphology and Luminosity 

The Central stellar cluster of the Milky Way has a stellar mass of approximately 4 x 1O 9 M0 [7] 
and harbors 2357 X-ray point sources within the central 17' x 17' alone [8]. In comparison, 
Omega Cen is the most massive stellar globular cluster in our Galaxy's halo with mass 
5 x 10 6 M Q [9], and has 45-70 X-ray sources associated with it [10]. 

The peak energy flux in gamma rays of the Central stellar cluster associated with the 
HG source is approximately 200 times that of Omega Cen, when corrected for distance. Given 
the three order of magnitude greater mass scale of the Central Galactic stellar cluster and its 
relative high-energy activity to that of the globular clusters, the factor of ~200 enhancement 
of flux is not unreasonable. 

In the same regard, the gamma ray profile of a globular cluster like Omega Cen only 
extends to 7.5', while that observed in the Central stellar cluster by HG is approximately 
30', which can also be accounted for by the 3 order of magnitude larger mass scale of the 
Central stellar cluster. 
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4 Conclusions 



The gamma-ray flux spectrum and profile observed by HG towards the Central stellar cluster 
of the Milky Way is consistent with other stellar cluster populations in observed globular 
clusters. The spectral index, exponential cutoff and peak flux energy is consistent with 
four of the eight detected globular clusters in the gamma-ray by Fermi-LAT. Therefore, 
there exists is no feature of the HG spectral source that necessarily requires a dark matter 
annihilation interpretation. 

Of the detected gamma-ray observed globular clusters, only 5 of 8 harbor detected 
MSPs [3]. The gamma-ray emission from stellar clusters such as the Central stellar cluster 
as well as globular clusters is expected to be dominated by MSP emission due to enhanced 
binary formation in these systems, where MSPs are spun up by accretion from their binary 
partners. Therefore, it is likely that the detection of the peaked spectrum towards the Galac- 
tic Central stellar cluster by Hooper & Goodenough is that of a population of millisecond 
pulsars bound within this massive stellar cluster. Obscuration, crowding, and absorption 
towards the central region of the Milky Way makes direct detection of the MSPs difficult, 
though this identification of the spectral signature of stellar cluster-like MSPs in the Galactic 
Center may motivate deeper searches. 
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